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CN , Abstract 

o 

Test particle geodesic motion is analysed in detail for the background spacetimes of the 
^-i degenerate Ferrari-Ibahez colliding gravitational wave solutions. Killing vectors have been 

used to reduce the equations of motion to a first order system of differential equations 
5JT), which have been integrated numerically. The associated constants of the motion have 

also been used to match the geodesies as they cross over the boundary between the single 
plane wave and interaction zones. 

X' 

U , 

1 Introduction 

Among the various theoretical arguments that make the study of gravitational waves one of the 
best arenas for probing classical and quantum gravity theories, those that involve the strongly 
nonlinear features of general relativity are very promising. The phenomenology related to the 
collision between two gravitational waves undoubtedly belongs to this category. However, 
even if the search for new exact solutions of the Einstein equations describing colliding waves 
essentially ended in the eighties (see e.g. |IJ for a review), the kinematics of test particle 
motion in such a background spacetime has never been studied in sufficient detail. 

In the literature it is known that the gravitational wave interaction generates a curva- 
ture singularity some time after the instant of collision and at a certain distance from the 
wavefronts. Moreover there exist degenerate solutions due to Ferrari and Ibahez [0] in some 
of which a horizon is created instead. The characteristics of geodesic motion in both types 
of degenerate solutions are studied in the present paper starting from some preliminary re- 
sults by Dorca and Verdaguer §. In particular we show that the horizon and singularity are 
reached in a finite time by those geodesies which fall into it. 

Two nontrivial Killing vectors are used to integrate the geodesic equations, reducing them 
to a first order system of differential equations which can then be integrated numerically. The 
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results are summarized in a series of appropriately chosen figures typical of the various kinds 
of behaviour that can occur. Finally some properties of the spacetime associated with the 
Riemann invariants, the Papapetrou fields and the principal null directions are discussed, em- 
phasizing the aspects arising from the asymmetry of these solutions relative to the coordinates 
spanning the surface of the wave front. 



2 The special class of Ferrari-Ibanez degenerate solutions 

In 1987 Ferrari and Ibahez H found two degenerate solutions of the Einstein equations that 
with an appropriate choice of the amplitude parameters associated with their strength can 
be interpreted as describing the collision of two linearly polarized gravitational plane waves 
propagating along z with opposite directions. Both these solutions are isometric to the interior 
of the Schwarzschild metric but one develops a singularity and the other a coordinate horizon: 

ds 2 = (1 + crsint) 2 (dz 2 - dt 2 ) + - — asmt ^x 2 + cos 2 z(l + crsint) 2 dy 2 , (1) 

1 + 0" sin t 

where a = ±1 (a = 1 for the metric with horizon and a = — 1 for the metric developing 
singularity). The interaction region where this form of the metric is valid (designated as 
"Region I" following 0]) is represented in the (t, z) diagram by an isosceles triangle whose 
vertex (representing the initial event of collision) can be identified with the origin of the 
coordinate system; the horizon/singularity is mapped onto the base of the filled isosceles 
triangle in Fig. 1. The only nonvanishing second order Riemann invariant for both cases is 
the Kretschmann scalar: 

48 

(1 +~o~smt) t 

In order to describe the larger spacetime of which this is only one region, one must 
introduce the two null coordinates 



Ki = „ | _.. 6 - (2) 



u = (t - z)/2 , v = (t + z)/2 t = u + v , z = v -u (3) 

in terms of which the metric (Q) takes the form 

ds 2 = — 4 [1 + o"sin(n + v)] 2 dudv 

l-asm(u + v) , 2 2 , , ri . . , l2 , 2 (4) 

H t r dx 2 + cos^ti - v)\l + <7sin(u + v) ] dy. 

1 + a sm(n + v) 

Following Khan-Penrose by proper use of the Heaviside step function H, one can easily 
extend the formula for the metric from the interaction region to the remaining parts of the 
spacetime representing the single wave zones and the flat spacetime zone before the waves 
arrive. The interaction region corresponds to the triangular region in the (u, v) plane bounded 
by the lines u = 0, v = and u + v = tt/2. One need only make the following substitutions 
in (1): 

u — > u H(u) v — ► v H(v) (5) 

that give rise to the four regions 



u > , v>0,u + v<ir/2 Region I Interaction region 

< u < tt/2 , v < Region II Single n-wave region 

u < , < v < tt/2 Region III Single v-w&ve region 

u < , v < Region IV Flat space 



(6) 
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Figure 1: The null coordinates and the different regions they induce. The (t, z) coordinates 
are scaled by a factor of \/2 with respect to the (u, v) . Region IV is flat space: geodesic paths 
therein are represented by straight lines whose slope affects the value of the constant K v in 
Region II. 



as shown in Fig. 1. In this way the extended metric in general is C° (but not C 1 ) along the 
null boundaries u = and v = 0. 

It is worth noting that any calculation concerning the spacetime metric can be appropriately 
done in u and v coordinates, while in the case of a test particle freely moving in that same 
background spacetime, calculations prove to be easier (especially in Region I) in the coordinate 
patch (t, z); so a switching between the two coordinate patches must always be kept in mind 
in reading what follows. 



3 Timelike geodesies 

For both degenerate solution background spacetimes we will solve the geodesic equations 

d 2 x a , , ,, dx^ dx v n ,_. 

^ + rV(*M)^^ = o (7) 

with the 4-velocity normalization condition specifying timelike orbits 

and where the coordinates, metric components and the Christoffel symbols depend on the 
proper time parameter r along the geodesies as explicitly indicated. Our goal here is to 
see how a massive particle freely transits from Region IV to Region II (or to its symmetric 
counterpart, Region III) and finally to Region I, as far as the chosen coordinate patches allow 
its geodesic motion to be described. Since the spacetime here admits 4 independent Killing 
vector fields, one can reduce the second order equations of motion to a first order system 
where three new constants (one for each Killing vector, two of them appearing always in a 
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Table 1: The Killing vectors in the different regions and the associated coordinates and con- 
stants they are referred to in the text. 



fixed analytic relation between each other: this leaves three) supplement the normalization 
condition above. One can then solve these four conditions for the first order derivatives 
= dx^/dr to obtain this system. 

For each Killing vector £ of the metric one obtains a constant of the motion from 

Z li U» = K (9) 

along each geodesic. In our case, for both kinds of metric, there are two obvious spacelike 
translational Killing vectors: £n\ = d x and £( 2 ) = 9 y , with associated constants K x = U ■ 
and K y = U ■ £(2)7 which always exist in both the single-wave and the interaction regions, 
since the metric (|l|) there does not depend on x or y. For the same reason, in Region II (or 
Region III) a third translational Killing vector is immediately found from the independence 
of the metric (||) on v (or on u), once expressed using (||); the associated constant will be 
designated by K v = U ■ £,S (or K u = U ■ It is also easy to see that two additional 

Killing vectors exist in Region I 

£(3) = cos y &z + siny tanz d y (10) 
£u\ = — siny d z + cosy tan z d y (11) 

Only one further constant is needed, which can be taken to be the square root of the sum of 
their two associated constants of motion 



Kz = \/K 2 z l+Kh (12) 

where K z \ = U ■ £/L and K z % = U ■ The timelike geodesies in Region I are then found to 
satisfy a subsystem of two first-order differential equations for motion in the t-z (or u-v) plane 
depending on three constants K x , K y , K z . For the sake of clarity, the Killing symmetries 
used are synthesized in Table 1. 

Note that by introducing the following differential operator: 

^ = (e(2)) 2 + (e { I 3)) 2 + (^4 ) ) 2 (13) 

and by letting it act on a function of the coordinates which oscillates along (x, y) with proper 
frequencies given by (K x , K y ), namely f(t, z, x, y) = e l ( Kx x+KyV ^ g(t, z), it is easy to see 
that, for any g(t, z), the following property is verified: 



£ 2 (f) = e 



i(K x x+K y y) 



^_ tallz Mt£)_4_ 9(t , £) 

OZ L oz cos z z 



(14) 
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which is clearly reminiscent of the square of the angular momentum operator, if one thinks at 
the specialization of the well-known Chandrasekhar-Xanthopoulos isometry || to the Ferrari- 
Ibahez degenerate solutions, that relates them to the "interior" Schwarzschild spacetime (i.e. 
for r <2m using standard Boyer-Lindquist coordinates). Actually it is exactly the operator 
associated with the integrable part of the Klein-Gordon equation (i.e. the z part), already 
studied by Dorca and Verdaguer and Yurtsever [Q, || : this gives the Killing vectors explicitly 
presented here a special importance and also helps to explain the integrability of the Klein- 
Gordon equation itself. 

To extend a geodesic from Region II to Region I the value of K z must be selected properly. 
The geodesic system in Region I is 



1 



dt{r) _ 
(It (1 + a sint(r)) 

dz(r) 1 



i Ki(l + g S mt(r)) K* + K* 



1 — (Tsint(r) (1 + (Tsint(r)) 2 



dr (1 + CTsint(T)) 2 
dx(r) 1 + a sint(r) 



'K* - Kl tan 2 z(r) 



dr 
dy{r) 



1 — a sin(r) 



1 



dr v cos 2 z(t) (1 + (Tsint(r)) 2 

while in Region II it is 



(15) 



du(r) 

dr 
dv(r) 

dr 
dx(r) 

dr 
dy{r) 

dr 



K v 1 
2 (1 +cr sin u(t)) 2 
1 



1 + r2 l + asmu(T) + r2 1 

x 1 — a sin u(t) y cos 2 u(t) (1 + a sin n(r)) 2 



K, 



2K V 
1 + 0" sin u(t) 
1 — a sin(r) 



1 



cos 2 u(t) (1 + cr smu(r)Y 



(16) 



Limiting the analysis to the (t, z) plane, we see that one of the most important differences 
between the horizon and singularity cases lies in the behaviour of the geodesies approaching 
t = it/2 (or u + v = it/2). 

• horizon case {a = 1) 

The slope of the trajectories in the limit t = it/2 is 

dz 



dt 







(17) 



(or du/dv = 1) for K x ^ 0, i.e. the test particle crosses the horizon with zero velocity 
in the direction of propagation of the waves, regardless for the values of K y and K z (see 
Fig. 5). Such behaviour is instead absent in the case K x = (see Figs. 3 and 7), where 
one finds 

dz 

" " (18) 



dt 



K*-K* tan 2 z 
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clearly depending on the values of K y and K z . The fact that making K x 7^ changes 
the behaviour of the geodesies at the horizon so much compared to the K x = case 
corresponds to a sort of "broken symmetry" between the x and y coordinates associated 
with the plane symmetry of these spacetimes, which we will briefly discuss in the next 
section. 

singularity case (a = —1) 

Either for K x = or K x ^ 0, we have: 



dz 



K 2 7 - Kl tan 2 z 

y (19) 



V K y + K z 



and the test particle approaches the singularity with z-velocity in general different from 
0. However, here it is important to distinguish whether K y vanishes or not. For K y = 
the geodesies exhibit a typical twofold behaviour (see Figs. 2 and 4) depending on 
whether they enter the interaction region or not. This is due to the existence of a 
critical value of K v , here denoted by -ft^ hr , in the single- wave region, or equivalently to 
a threshold for the velocity along the z axis in the vacuum region. Geodesies in Region 
II with K v = K^ r can be extended into Region I with K z = 0, i.e. they remain at 
a fixed z position, irrespective of the presence of the waves; this K^ hr and the Region 
II-Region I crossing value u* of u are related by: 



K^ = (sin u* - 1) Jl + Kl . (20) 

V 1 + sin u* 

For particles exceeding this value of K v , the geodesies are of Type I, i.e. entering the 
interaction region and reaching the singularity with dz/dt — > 1, otherwise, they are of 
Type II, i.e. confined in the single-wave region and approaching the singularity with 
dz/dt -> -1. 

As stated above, along the geodesic corresponding to the critical value of K v , the particle 
enters at a certain position the interaction region with dz/dt = and doesn't move. 
This geodesic marks a triangular part of Region I, also delimited by the singularity and 
the u axis, which actually doesn't allow any geodesic trajectory to enter. 

For K y ^ the threshold disappears, as shown in Fig. 6. 

As for the proper times of those geodesies entering the interaction region and approaching 
the line t = vr/2, either in the horizon or singularity cases, the situation is summarized in 
Tables 1 and 2, for the case of a particle entering Region II at v = —vr/2, with arbitrary 
choices of K x , K y and giving K v the value that corresponds to zero initial velocity along z 

(i.e., du/dv = 1 for u — > 0), which is — J K% + K | + 1 . 

4 Papapetrou fields and principal null directions in Region I 

The Killing vector field generates a Papapetrou field F whose principal null directions 
are aligned with those (repeated) of the spacetime itself (which is of Petrov type D). The key 
properties are the following: 

E(X) afl = C af3fMV X? X v = -3 (1 - a sin t) \a(X) ® a(X)}^ , 

H{X) atl = *C aPlxv XPX v = Q (21) 
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Table 2: (Horizon case) the proper times for different initial data sets. 
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Table 3: (Singularity case) the proper times for different initial data sets. (-): the geodesic 
doesn't reach Region I. It is evident that, while the particles approach the horizon smoothly, 
the formation of a singularity makes their transit in Region I almost immediate. 



where X = £(i)/||£(i)|| = ^/gl^dx, and 

a{X) = V x X = -[d lx \n x /g x ^]dx» = — dt (22) 

COS % 

is the (timelike) acceleration field associated to X and (TF) stands for the trace free part of 
a tensor. Explicitly, for a generic Q) tensor A, orthogonal to X, one has: 

A aP = A «P ~ \ \^ A \ P ( X U . p (X) =g-X®X, (23) 

with P(X) being the projector orthogonal to the spacelike X direction. Finally the Papa- 
petrou field associated to £ri) is 

F a /3 = Z(l)a;f> = K(l) A a(X)] af3 . (24) 
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The two independent principal null directions of both F a p and C a p^s are 



I = 



1 



2 cost 
cost 



(1 + asint) 2 
Ot H : Or 



n 



(1 + fjsint) 2 



cost 
d t - d x , 



(25) 



affinely parametrized and normalized so that I ■ n = — 1. From these expressions it is clear 
that the special role of the Killing vector (compared to £(2)) is due to the fact that it 
belongs to the 2-plane spanned by the principal null directions I and n. This, at least from 
a purely geometrical point of view, leads to a difference between the x and y coordinates in 
terms of which the metric is expressed. 



5 Concluding remarks 

A detailed analysis of (timelike) geodesic motion has been carried out in the spacetimes of 
the Ferrari-Ibahez degenerate solutions representing two colliding gravitational plane waves. 
The study of these simple cases is a useful exercise in general relativity because they contain 
the main features related to the nonlinear collision of waves, i.e., the creation of a horizon or 
that of a curvature singularity. 
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Figure 3: For a = 1 (metric with horizon) and = 0, K y = 0), the geodesies regularly fill 
the entire Region I, coming from a set of initial data characterized by a limiting value of K v , 
K^ m = —1.6653 under which the particles possess a too large momentum along z (moving 
into reverse with respect to the single wave) to enter the collision region. 
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Figure 5: For a = 1 and {K x = 1, K y = 0), the geodesic motion has a global homogeneity for 
all values of K v (it is always the case when K x ^ 0). As one expects, in Region II the particles 
are accelerated (in the sense that they exhibit a relative acceleration d 2 z/dt 2 ) in the positive 
z direction, while in Region I they are decelerated until they reach zero velocity exactly at the 
horizon. The limit K v — > — oo reconstitutes the so-called fold singularities already observed 
in the case of null geodesies by Dorca and Verdaguer Q 
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Figure 6: The case a = — 1 and (K x = 3, K y = 1) is a good example of what the geodesies of 
the metric with singularity look like when K y ^ 0. The fold singularities arise for K v — > — oo 
while the particles accelerate abruptly (relative acceleration of course) and with a change in 
sign along z in Region II. There are no forbidden regions in Region I. 
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Figure 7: For a = 1 and (i^ = 0, = 1) (this case qualitatively characterizes all the sets 
of initial data with K x = 0), the metric with horizon exhibits a very homogeneous geodesic 
motion in Region II, where all the possible values of the momentum along z are allowed to 
enter Region I (so generating again the fold singularities) and slight changes in Region I, 
where the particles decelerate for high entering values of dz/dt and accelerate otherwise. 
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